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Abstract 

We evaluate the gravitino production rate in supersymmetric QCD at high tem- 
perature to leading order in the gauge coupling. The result, which is obtained by 
using the resummed gluon propagator, depends logarithmically on the gluon plasma 
mass. As a byproduct, a new result for the axion production rate in a QED plasma 
is obtained. The implicatons for the cosmological dark matter problem are briefly 
discussed, in particular the intriguing possibility that gravitinos are the dominant 
part of cold dark matter. 



1 Introduction 



Supersymmetric theories, which contain the standard model of particle physics and 
gravity, predict the existence of the gravitino p], a spin-| particle which aquires a mass 
from the spontaneous breaking of supersymmetry. Since the couplings of the gravitino 
with ordinary matter are strongly constrained by local supersymmetry, processes involving 
gravitinos allow stringent tests of the theory. 

It was realized long ago that standard cosmology requires gravitinos to be either very 
light, ttlq < 1 keV [2], or very heavy, ttlq > 10 TeV [3J. These constraints are relaxed if the 
standard cosmology is extended to include an inflationary phase [HE]. The cosmo logically 
relevant gravitino abundance is then created in the reheating phase after inflation in which 
a reheating temperature Tr is reached. Gravitinos are dominantly produced by inelastic 
2^2 scattering processes of particles from the thermal bath. The gravitino abundance 
is essentially linear in the reheating temperature Tr. 

The gravitino production rate depends on m^/m^, the ratio of gluino and gravitino 
masses. The ten 2^2 gravitino production processes were considered in [5J for m g <C itiq. 
The case m g 3> tuq, where the goldstino contribution dominates, was considered in [6]. 
Four of the ten production processes are logarithmically singular due to the exchange of 
massless gluons. As a first step this singularity can be regularized by introducing either a 
gluon mass or an angular cutoff |5J. The complete result for the logarithmically singular 
part of the production rate was obtained in [7]. The finite part depends on the cutoff 
procedure. 

To leading order in the gauge coupling the correct finite result for the gravitino pro- 
duction rate can be obtained by means of a hard thermal loop resummation. This has 
been shown by Braaten and Yuan in the case of axion production in a QED plasma 
[8]. The production rate is defined by means of the imaginary part of the thermal axion 
self-energy [9]. The different contributions are split into parts with soft and hard loop 
momenta by means of a momentum cutoff. For the soft part a resummed photon propaga- 
tor is used, and the logarithmic singularity, which appears at leading order, is regularized 
by the plasma mass of the photon. The hard part is obtained by computing the 2^2 
scattering processes with momentum cutoff. In the sum of both contributions the cutoff 
dependence cancels and the finite part of the production rate remains. For the gravitino 
production rate the soft part has been considered in [10J and the expected logarithm of 
the gluon plasma mass has been obtained. 

Constraints from primordial nucleosynthesis imply an upper bound on the gravitino 
number density which subsequently yields an upper bound on the allowed reheating tem- 
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perature Tr after inflation [H]-[l5]. Typical values for Tr range from 10 7 — 10 10 GeV, 
although considerably larger temperatures are acceptable in some cases [16]. In models of 
baryogenesis where the cosmo logical baryon asymmetry is generated in heavy Major ana 
neutrino decays [T7j, temperatures Tr ~ 10 8 — 10 10 GeV are of particular interest |18j . 
Further, it is intriguing that for such temperatures gravitinos with mass of the electroweak 
scale, i.e. uIq ~ 100 GeV can be the dominant component of cold dark matter [7J. In 
all these considerations the thermal gravitino production rate plays a crucial role. In this 
paper we therefore calculate this rate to leading order in the gauge coupling, extending a 
previous result [7J and following the procedure of Braaten and Yuan [8] . 

The paper is organized as follows. In section 2 we summarize some properties of 
gravitinos and their interactions which are needed in the following. In order to illustrate 
how the hard thermal loop resummation is incorporated we first discuss the axion case 
in section 3. The most important intermediate steps and the final result for the gravitino 
production rate are given in section 4. Using the new results the discussion in [7J on 
gravitinos as cold dark matter is updated in section 5, which is followed by an outlook in 
section 6. The calculation of the hard momentum contribution to the production rates is 
technically rather involved. We therefore give the relevant details in the appendices. 

2 Gravitino interactions 

In the following we briefly summarize some properties of gravitinos which we shall 
need in the following sections. More detailed discussions and references can be found in 

PU EDI EHj. 

Gravitinos are spin-3/2 particles whose properties are given by the lagrangian for the 
vector-spinor field ijj^(x), 

£ = -\e^»lslvd P ^ ~ \m^,[Y, 7*1^ - ■ C 1 ) 

Here vtiq is the gravitino mass, M = {ZttGn)- 1 ' 2 is the Planck mass and S M is the 
supercurrent corresponding to supersymmetry transformations. W an d are Majorana 
fields, so that = S^. 

Free gravitinos satisfy the Rarita-Schwinger equation, 

- V p<T 757^ CT - \m Q [f, 7>, = , (2) 

which, using 

= , &M x ) = > ( 3 ) 

3 



reduces to the Dirac equation 



(ip-me)il> lt (x) = 0. (4) 

Consider as matter sector first a non-abelian supersymmetric gauge theory with la- 
grangian 

C = -i J£,JW + {vfD^ - m- g 5 ab ) \ b (5) 

for the vector boson A ^ and the gluino A a . Supersymmetry is explicitly broken by the 
gluino mass term. Hence, the supercurrent is not conserved, 

= m-gS . (6) 

The calculation of the gravitino production rate in section 4 will involve squared matrix 
elements which are summed over all four gravitino polarizations. The corresponding 
polarization tensor for a gravitino with momentum P reads 

I 

= -<? + m ,) (,„ - p -^)-\ (v + %)<r- m 6) (r + . <T) 

Since ip^x) is a solution of the Rarita-Schwinger equation one has for the polarization 
tensor 

YU^(P) = 0, PTV(P) = 0, (8) 

(P-m)U^(P) = 0. (9) 

We shall be interested in the production of gravitinos at energies much larger than the 
gravitino mass. In this case the polarization tensor simplifies to 

iV(p) ~ -f 9lw + h^f , (io) 

where we have used 7^5^ = 0. Clearly, the first term corresponds to the sum over the 
helicity ±| states and the second term represents the sum over the helicity ±~ states 
which represent the goldstino part of the gravitino. Eq. (TTU|) is the basis for the famil- 
iar substitution rule ip^ — > yf^r^tV' which is used to obtain the effective lagrangian 
describing the interaction of goldstinos with matter 
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The gravitino production rate at finite temperature can be expressed in terms of the 
imaginary part of the gravitino self-energy [9] which takes the form, 



E(P) 



tr 



lWP)£g*(P) 



n Mi/ (p)^(p)...^(P) 
(-p)^(P)...^(P) 



2m? 



tr 



PS(P) . . . S(P) 



(11) 



Here we have used eq. ([6]) for the divergence of the supercurrent. Note, that F = ^/3m^M 
is the scale of spontaneous supersymmetry breaking, which gives the strength of the 
goldstino coupling to the supercurrent. The dots denote the sum over the contributions to 
the self-energy in the loop expansion. In S^(P) = j[y u , , j p }j' J '\ a (Pi)F up (K 1 ) and S(P) = 
j [7", 7 p ]A a (Pi) F vp (Kx) , with P = P\+K\, A a (Pi) represents one end of an internal gluino 
line and F yp (Ki) stands for the end of one or two internal gluon lines. For uiq <C m- g the 
goldstino part dominates the gravitino production cross section. 

Eq. (fTTj) is useful to derive relations between the helicity ±| and the helicity ±| 
contibutions to the self-energy. As shown in Appendix A, one obtains to two-loop order 



S(P) oc 



M 2 



1 + 



m 9 
3 ™L 



(12) 



We shall exploit this fact to perform the hard thermal loop resummation, which is neces- 
sary because of the infrared divergences, just for the helicity ±| part of the rate. 

The full supercurrent also involves quarks and squarks in addition to gluons and 
gluinos. The divergence of the additional part of the supercurrent that involves the cubic 
gravitino-quark-squark coupling is again proportional to the parameter of supersymmetry 
breaking, i.e. m~, the squark mass squared. The corresponding goldstino contribution 
to the production rate is then proportional to m|, which is suppressed at high energies 
compared to the gluino contribution for dimensional reasons. 



3 Axion production 

Let us now consider the thermal production of axions in a relativistic QED plasma of 
electrons and photons. According to the procedure outlined in the introduction the ther- 
mal production rate can be obtained as sum of two terms, a soft momentum contribution 
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which is extracted from the axion self-energy evaluated with a resummed photon propa- 
gator and a hard momentum contribution which is computed from the 2 — > 2 processes. 



Y 




Figure 1: Axion self energy; the blob denotes the resummed photon propagator. 



The axion-photon interaction is described by the effective lagrangian 

C = -Jf aF ^ plJV ' ( 13 ) 

where / is the axion decay constant. The axion self-energy U a (P) (cf. fig. [1]) can be 
evaluated in the imaginary-time formalism for external momentum P = (po>p)> with 
Po = i2imT and p = |p|. In covariant gauge the resummed photon propagator has the 
form [251 EE] 

iA^(K) = % (A^At + B^A L + C^) , (14) 



with the tensors 



1 



-V = ~9iiv~ 7a K v^- K •v(K tl v v + KvV^+K^Ky 



k 2 

K-v /K-v^ 2 



K K u 

and the transverse and longitudinal propagators 
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A T (k ,k) 



k 2 -k 2 -U T (k ,k) ' 



A L (k ,k) = k2 J L{k ^ ky (16) 

Here K = (fco,k) with ko = i2imT and k — |k|, £ is a gauge-fixing parameter, v is the 
velocity of the thermal bath and Ht/l are the transverse and longitudinal self-energies of 
the photon. The corresponding propagators At/l have the spectral representation 

r°° 1 

A T/L (k ,k)= du- Pl/t(u, k) . (17) 

J — oo Kq — 00 



For \u\ < k the spectral densities Pl/t are given by [26] . 
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p T (aJ,k) 



X 



Am 2 (1 - x 2 )(A T (x) 2 + (z + B T (x)) 2 ) 
3 2x 



Am 2 A L (x) 2 + (z + B L (x)) 



2 ; 



(18) 



where m 1 = eT/3 is the plasmon mass of the photon, x = u/k, z = k 2 jmt and 



3 3 / x 2 

At{x) = -nx, B T (x) = - 2- + x In 

4 4 y 1 — ar 

A L (x) = ^nx, B L {x) = ^( 2_xln 3~~ 



1 + x 
1 — X 



(19) 



The contribution to the axion production rate from soft virtual photons is obtained by 
analytically continuing the axion self energy function U a (P) from the discrete imaginary 
value p to the continuous real value E = p [8] , 



Imn a (E + ie,p) 



^J" 1 /*^cut 



87T/ Z JO J-k UJ 



k duo 



Pl(w, k) (l-^)+ Pt(uj, k) M - ^ J 



.(20) 



The axion production rate depends logarithmically on k cut . The corresponding coefficient 
can be obtained analytically. The remaining constant has to be evaluated numerically. 
This yields the result, first obtained in [8], 



ft/ 3m 2 T 

K oit (E) - 



ln^b - 1.379 



mt 



16tt/ 2 

The corresponding collision term in the Boltzmann equation is 

d 3 p 



(21) 



c s a oit (T) = J 



where 



(2tt) 3 
e 2 C(3)T 6 
24vr 3 / 2 

n B (E) 



n B {E)Tf\E) 

fx r 



In 



^cut 



- 0.689 



exp (E/T) - 1 



(22) 



(23) 



is the Bose-Einstein distribution. 



7 




Figure 2: Axion production in electron-photon scattering. 



The dependence of r^ oft (P) on the cutoff k cut is cancelled by the contribution from 
hard virtual photons to the self-energy. This part of the axion production rate can be 
obtained directly from the processes e ± 7 — > e ± a (cf. fig. [2]) [9], 



n B (E)T^ A (E) 



dVt p 1 
17 2P 



n 



d 3 Pi 



[ (27r) 3 2^_ 



(2tt) 4 5 4 (P 1 + P 2 -P-P 3 



n F (E 1 )n B (E 2 )(l - n F (E 3 ))\M\ 2 &(\ Pl - p 3 | - k, 



cut J 



(24) 



Here we have used rotational invariance by averaging over the directions of the axion mo- 
mentum; |M| 2 is the photon-axion matrix element squared for e~(Pi)7(P2) —* ^ (Ps)cl(P), 



\M\ 2 = f 2 [-^-2s-t 



2s_ 2 
t 



(25) 



where s = (Pi + P2) 2 , t = (Pi — P3) 2 , and n F (E) is the Fermi-Dirac distribution, 



n F (E) 



exp(P/T) + 1 



(26) 



The phase space integration has to be carried out under the constraint on the virtual 
photon momentum A; = |pi — P3 1 > k cnt . For the angular integrations it turns out to be 
convenient to define all momenta with respect to k. Some details of this calculation are 
given in appendix B. One finally obtains, 

3e 2 1 



n B (E)TT & (E) 



2 8 7T 3 / 2 P 2 

X 



dEidE 3 dkn F (Ei)n B (E 2 )(l - n F (E 3 )) 



(Ei - P 3 ) 2 - k 2 



(E 3 + Ei) 2 (E + E 2 ) 
k* 



, 2 E\ + E\ + 2EE 2 
T H 

3 k 2 



n 



(27) 



where E2 = E + E 3 — E\ and the integrations are restricted by Jl, 

Q = &(k - k cut )Q(k - \E, - E 3 \) 

Q(E 1 + E 3 - k)Q(2E + E 3 -E 1 - k) 

e(E 1 )e(E 3 )e(E + e 3 - e x ) . 
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After performing the /c-integration one is left with several domains for the E\- and £3- 
integrations. The logarithmic dependence on k CVLt can be extracted by means of a partial 
integration in E\. In the remaining part of the integral k cu t can be set equal to zero. The 
final result reads 



-ihard 



(E) 



2vr 2 



-T 6 In 



2T 

k rn t 



17 

t~ 7 



167T 3 / 2 I 

+ (e E ' T - l) jf dE 3 [1 - n F (E 3 )] j Q dE 
\E\ — E 3 



cm) 

E+E 3 

1 1: 



x In 



E 3 



x 



Q(E - E 1 



S 

E 2 



n F (E 1 )n B (E 2 



2E X 



2 ~t + \ Et 



E 3 E 2 
El 



-Q{E X -E 3 )^ (2E! 



E- 



^ + (El 



E 



n F {E{) + n B (E 2 ) \ 
T ) 
n F {E 1 ) + n B (E 2 ) 



T 



+Q(E 3 -E 1 )[2E 1 + (E 2 1 +El 



n F {E x ) + n B {E 2 ) 



(29) 



This result agrees with the one obtained in [Hj except for the first expression cx Q(E — E{) 
in the double integral. Integration over the axion energy E yields for the collision term 
in the Boltzmann equation 

d 3 p 



C hard (T) = J 



(2tt) 3 
e 2 C(3)T 6 
24vr 3 / 2 



n B {E)T^\E) 
2T 



In 



17 _ ('(2) 
+ 6 7+ C(2) 



1.280 



(30) 



The numerical constant is about 20% smaller than the one obtained from the axion rate 
given in [8]. 

Consistency requires that the dependence on the cutoff k cnt cancels in the total pro- 
duction rate. Comparison of eqs. fl22l and fT29|) shows that this is indeed the case. The 
result for the total axion collision term reads 



C a (T) = C s a ott (T) + C™ Id (T) 



2 C(3)T 6 



48vr 3 / 2 



In 



2 



0.8194 



7- 



(31) 
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4 Gravitino Production 



The rate for the thermal production of gravitinos can be calculated in complete analogy 
to the axion production rate. It is dominated by QCD processes since the strong coupling 
is considerably larger than the electroweak couplings. The contribution due to soft virtual 
gluons can again be extracted from the gravitino self-energy with a resummed gluon 
propagator to which the contribution from hard 2 — > 2 processes has to be added. 

Properties and interactions of the gravitino have been discussed in section 2. For 
supersymmetric QCD with gluons, gluinos, quarks and squarks one obtains [TH [20j [21] , 



i 



V2M 



(D^^YYPlX - (D^)xPr1 v 1^v\ - m ^ \l\ 71 1^K P - (32) 



i 



Here \ denotes a left-handed quark or antiquark and the corresponding squark. For 
light gravitinos one can use a simpler effective lagrangian [221 1231 [21] . The corresponding 
goldstino-gluon-gluino coupling can be read off from eqs. (jSJ) and ( TlTT) 



777 ~ 

C cS = - « V [7", 71 ^ a FL + • • • • (33) 



2^6 Mm d 

Here ip is the goldstino, the spin- 1/2 component of the gravitino. The effective theory 
contains the same vertices as the full theory, except for the gravitino-quark-squark-gluon 
vertex. Instead, there is a new four particle vertex, the gravitino-gluino-squark-squark 
vertex [21]. All vertices are proportional to supersymmetry breaking mass terms, i.e., m~ 
and trig. At high energies and temperatures, with rrig, m g -C T, contributions involving 
the cubic goldstino-quark-squark coupling are suppressed by m?/T 2 relative to the gluino 
contribution because of the higher mass dimension of the coupling. 




Figure 3: Gluon-gluino loop diagram, the leading contribution to the imaginary part of 
the gravitino self energy. The blob denotes a resummed gluon propagator. 
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In the imaginary-time formalism one obtains for the goldstino self-energy (cf. fig. [3]) 
with momentum P summed over helicities: 



E(P) = tr 



l=±l/2 



tJ2 

k 



d 3 k 



rtr 



PV»-A^(K)V{ 



(34) 



Here we have neglected gluino and gravitino masses since rriq, m g -C T; V^ 2 are the 
vertices, Q = P — K is the momentum of the gluino, and A^^K) is the resummed 
gluon propagator, which is obtained from the resummed photon propagator (fT4"l) by the 
substitution m 7 — > m g . The thermal gluon mass for N colours and rif colour triplet and 
anti-triplet chiral multiplets is given by 



9 2 T 2 



<N + n f ). 



(35) 



9 Q 

This result is easily obtained from the expressions for the gluon vacuum polarization [2 
by adding up the contributions from gluons, gluinos, quarks and squarks. 

Inserting gluon propagator and vertices in eq. (I3"4"l) yields the gauge-independent result 

d 3 k 



4 m 2 -T 

G k 



(27T) 



(D L Aj 



D T A T ) - 



(36) 



where 



D T (k ,k,E,p, pk) 



1 

32 
1 



tr{f $,Y}®W,^}A^} 



D L (k ,k,E,p,pk) = -tr{f $,Y]@$,^]B^} . 



(37) 



After a straightforward calculation, analogous to the one for the axion self-energy, one 
finds for the gravitino production rate 

ImS(P + ie,p) 



soft 

1 G 



E 

m 2 g{N 2 - 1)T 
6nM 2 m% 

Gr 



h <[/ccut 

^cut 







k UJ 



UJ 



Pl(w, k) 1 - — + p T (uJ, k) 1 - — 



A' 2 



(3? 



The momentum integral depends logarithmically on the cutoff k CVLt . The integrand, which 
is identical with the one for the axion production rate (1201) . agrees with the result obtained 
in [10] . After performing the momentum integrations one finally obtains 

(N 2 - \)m\m 2 g T 



TH Soft 

1 G 



(E) 



4nM 2 m% 

Gr 



1.379 



(39) 
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Note, that the overall normalization differs from the expression given in [10] by the factor 
4(iV 2 - 1). 

The dependence of the soft part of the gravitino production rate on the cutoff fc cu t 
is again cancelled by the cutoff dependence of the contribution from the hard 2^2 
processes. There are 10 processes denoted by A to J [5]: 



A: g a + g b -> g c + G 



,G g a 



+ 



b ,<F 



, ~c 
J g 



g 




G 



+ 




G g a 




G 



~c 

g 



B: g a + g b — ► g c + G (crossing of A) 
C: qi + g a ^ q 3 + G 



G q, 



1.i 




G % 
+ 

a 

q § 




G q, N G 
+ 

1j g " X M; 




D: g a + qi — > 9j + G (crossing of C) 
E: ^ + qj — > g a + G (crossing of C) 
F: g a + g b ^ g c + G 




g 



+ 



g g 




G: q, + g a -> ^ + G 




+ 



G 



+ 




H: qi + g a ^ q d + (5 
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• I: qi + qj — > g a + G (crossing of G) 

• J: qi + qj — > g a + G (crossing of H) 

The corresponding matrix elements have been evaluated in [7j. As discussed in section 2, 
they must have the form 

^ « w f 1 + 4) (40) 

in the high energy limit. 





process i 


\Mi?/& (1+4! 


) 


A 


g a + g b -^g c + G 


A(s + 2t + 2^)\f abc \ 


2 


B 


g a + g b ^g c + G 


-A(t + 2s + 2^)\f al 


c 2 


C 


q~i+g a -> qj + G 




D 


g a + q t -> qj + G 


-2t|^| 2 


E 


Qi + qj^g a + G 


-2t\T^ 


F 


g a + g b ^g c + G 




} 2 +st+t 2 ) 2 I rabc 12 
st(s+t) \J 1 




G 


qi + g a -+ qj + G 


-4(, + ^)|^| 2 


H 


q~i + g a -> Qj + G 


-2(t + 2s + 2^)\T- 


2 


I 


Qi + Qj -> g a + G 


-4(t+?)|^| 2 


J 


Qi + Qj -> g a + G 


2(s + 2t + 2f)|T^| 2 





Table 1: Squared matrix elements for gravitino (G) production in two-body processes in- 
volving left-handed quarks (qi), squarks (qi), gluons (g a ) and gluinos (g a ). The values are 
given for the specified choice of colors and summed over spins in the initial and final state, 
jabc an ^ rpa^ are ^ g usua i SU(3) colour matrices. 
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In table [T] the squared matrix elements of all ten processes are listed. Sums over initial 
and final spins have been performed. For quarks and squarks the contribution of a single 
chirality is given. One easily checks that the matrix elements satisfy the relevant crossing 
symmetries. The particle momenta P\, P 2 , P3, and P used in the calculations correspond 
to the particles in the order in which they are written down in the column "process z" of 
table [H This fixes the energies of Bose and Fermi distribution. The matrix elements in 
the table correspond to the definitions s = (Pi + P2) 2 and t = (Pi — -P3) 2 . 

The different processes fall into three classes depending on the number of bosons and 
fermions in initial and final state. A, C and J are BBF processes with two bosons in the 
initial and a fermion in the final state; correspondingly, B, D, E and H are BFB processes, 
and F, G and I are FFF processes. Only four processes, B, F, G and H contribute to the 
logarithmic cutoff dependence. The gravitino production rate is then given by (cf. [9J), 



n F (E)T h f d (E) 



dVtp 1 
17 2E 



-pr d 3 pj 

L\ ^f2E t 



(2tt) 4 5 4 (P 1 + P 2 - P - P 3 



x [n BBF \M BBF \ 2 + n BFB I M BFB | 2 + n FFF \M FFF \ 2 ^ 6(|pi - p 3 | - fc cut ) . (41) 

Here, n BBF , n BFB and n FFF are the products of number densities for the corresponding 
processes, e.g., 

n BBF = n B (Ei)n B (E 2 )(l - n F (E 3 )) . (42) 

The matrix elements \M BBF \ 2 etc. are obtained by summing the corresponding matrix 
elements in table [1] with the appropriate multiplicities and statistical factors. Angular and 
momentum integrations can now be carried out as in the case of axion production. One 
finally obtains the result 



■phard 



(E) 



1 



3m% 



8tt 3 M 2 
+ (N + n f ) (e E ' T + l 



hut J 6 1 C(2), 



x 



Q(E - Ei) 
-Q(Ei - E 3 
+&(E 3 - Ei 



d 



dEi 

, d 



00 f E+E 3 f\E _ ft 

dE 3 dEilnl 1 1 3 

n JO \ E 3 

7?2 77i2 

(n BFB + n FFF )( 



dEi 
d 



E 2 

(n BFB + n FFF )-^(El + El) 



dEi 

+Ibbf + Ibfb + Ifff ( ■ 



(n BFB + n FFF ) (E\ + El 



(43) 
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Performing the differentiations with respect to E\ yields expressions analogous to the one 
given in eq. (1291) . Ibbf, Ibfb and Ifff, which are not all proportional to N + Uf, are 
given in appendix C; they contribute to the cutoff-independent part of r^ ard . 

The dependence on k cut cancels in the sum of soft and hard contributions to the 
production rate. From eqs. fl39|) and f|43|) one obtains for the collision term 

C d {T) = I ^n F (E)(r^E) + T h - d (E) 




(N + rif) + 0.5781^/ \. (44) 



This is the main result of this paper. It allows to calculate the gravitino abundance to 
leading order in the gauge coupling g(T), contrary to previous estimates which depended 
either on ad hoc cutoffs [HI [6] or on an unknown scale of the logarithmic term [7] . 

An important question concerns the size of higher-order corrections. Note, that g ~ 
0.85 for T ~ 10 10 GeV. This is much better than at the electroweak scale T ~ 100 GeV 
where g ~ 1.2, or for the quark-gluon plasma at T ~ 1 GeV where g ~ 2.5. However, 
one still has to worry about the usually assumed separation of scales g 2 T <C gT <C T, 
which would correspond to fi g <C m g <C T, where \i g ~ g 2 T is the magnetic screening 
mass. Note, that for the supersymmetric standard model with N = 3 and rif = 6 one has 
m g ~ T. For the static Debye and magnetic screening masses the separation of scales has 
recently been studied in detail for the case of non-supersymmetric QCD [28]. For real- 
time processes almost nothing is presently known about non-perturbative effects related 
to the magnetic screening mass. This is a challenging theoretical problem. 



5 Gravitinos as cold dark matter 

We can now study the cosmological implications of our result eq. f|44|) for the Boltz- 
mann collision term of gravitino production. We are particularly interested in the case of 
large reheating temperatures after inflation, i.e. Tr ~ 10 8 — 10 10 GeV, which are relevant 
for models of leptogenesis. In the following we shall concentrate on the possibility that 
the gravitino is the lightest supersymmetric particle (LSP), updating the discussion in [7], 
where it was pointed out that a large gravitino mass ~ 100 GeV is compatible with 
such reheating temperatures. We shall ignore the non-thermal production of gravitinos 
[29| [30] which depends on the model of inflation. 
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From the Boltzmann equation 

dn 



+ 3Hn d = C 6 , (45) 

one obtains for the gravitino abundance at temperatures T < Tr, assuming constant 
entropy, 

Y (T)= n ^ T) ~ 9 ^ T) °^ Tr) (46) 
G{ ) n vad (T) 9,s(Tr) H{T R )n vaA {T R ) ' 1 ' 

where g*s(T) is the number of effectively massless degrees of freedom [21]. For T < 1 MeV, 
i.e. after nucleosynthesis, g+s(T) = ||, whereas g±s(T R ) = 215 j n the supersymmetric stan- 
dard model. With H(T) = (g*(T)ii 2 /90y^ 2 T 2 /M one obtains in the case of light gravitinos 
(tjiq -C rrig(fi), p ~ 100 GeV) from eqs. (j4"6j) and (jUj) for the gravitino abundance and 
for the contribution to Qh 2 , 

Y — 11- itr" ( — 1 f 10 ° GeV ^ 2 2 (47) 

fi^/i 2 = m (5 F (5 (T)n rad (T)/i 2 p c r 1 

" °- 2 Hib^G^J(,^^J(,TiwJ • (48) 

Here we have used g(T R ) = 0.85, n ra d(T) = ((3)T 3 /tt 2 , and m g {T) = g 2 (T) / g 2 (p)m g (p)] 
p c = 3HqM 2 = 1.05/i 2 10~ 5 GeV cm -3 is the critical energy density. The new result for 
VLqK 2 is smaller by a factor of 3 compared to the result given in [7J. Due to the large value 
of the plasma mass m g an estimate of the gravitino production rate, which is based just 
on the logarithmic term of the 2 — > 2 cross sections as in [7J is rather uncertain. 

It is remarkable that reheating temperatures Tr ~ 10 8 — 10 10 GeV lead to values 
Qgh 2 = 0.01 ... 1 in an interesting gravitino mass range. This is illustrated in fig. H]for 
a gluino mass m g = 700 GeV. As an example, for Tr ~ 10 10 GeV, tuq ~ 80 GeV and 
h ~ 0.65 [21] one finds Qq = 0.35, which agrees with recent measurements of Qm |31j . 
In general, to find a viable cosmological scenario one has to avoid two types of grav- 
itino problems: For unstable gravitinos their decay products must not alter the observed 
abundances of light elements in the universe, which is referred to as the big bang nucle- 
osynthesis (BBN) constraint. For stable gravitinos this condition has to be met by other 
super particles, in particular the next-to-lightest super particle (NSP), which decay into 
gravitinos; further, the contribution of gravitinos to the energy density of the universe 
must not exceed the closure limit, i.e. Qq = Pq/ p c < 1- Consider first the constraint 
from the closure limit. The condition Qq = YqItiqUt-^/ p c < 1 yields an allowed region 
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Figure 4: The density parameter Q^h 2 for different gravitino masses itlq as function of 
the reheating temperature Tr. The gluino mass has been set to m g = 700 GeV. 



in the m^-m^ plane which is shown in fig. [5] for three different values of the reheating 
temperature Tr. The allowed regions are below the solid lines, respectively. 

With respect to the BBN constraint, consider a nonrelativistic particle X decaying 
into electromagnetically and strongly interacting relativistic particles with a lifetime tx- 
X decays change the abundances of light elements the more the longer the lifetime Tx 
and the higher the energy density mxYxn Ta d are - These constraints have been studied in 
detail by several groups [HI [121 HH]- They rule out the possibility of unstable gravitinos 
with m 6 ~ 100 GeV for T R ~ 10 10 GeV. 

For stable gravitinos the NSP plays the role of the particle X. The lifetime of a fermion 
decaying into its scalar partner and a gravitino is 



TNSP 



m%M 2 
48?r^ 



in 



(49) 



NSP 



For a sufficiently short lifetime, r^sp < 2 • 10 6 s, the energy density which becomes free in 
NSP decays is bounded by m x Yx < 4 • 10~ 10 GeV, which corresponds to flxh 2 < 0.008. 
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Figure 5: Upper and lower bounds on the gluino mass and the NSP mass as functions of the 
gravitino mass. The full lines represent the upper bound on the gluino mass m- g > m^sp 
for different reheating temperatures from the closure limit constraint. The dashed line is 
the lower bound on m^sp which follows from the NSP lifetime. 



The lifetime constraint yields a lower bound on super particle masses which is represented 
by the dashed line in the mg-mNSP/g plane in Fig. [5j 

In order to decide whether the second part of the BBN constraint, f^NSP^ 2 < 0.008, 
is satisfied, one has to specify which particle is the NSP. The case of a higgsino-like 
neutralino as NSP has been discussed in [7]. A detailed discussion of the case where a 
scalar r-lepton is the NSP has been given in [T3].|15|. 

A complete treatment of gravitinos as cold dark matter has to include non-thermal 
contributions. The situation is analogous to leptogenesis where, in principle, non-thermal 
contributions also have to be added to the thermal part. However, non-thermal contribu- 
tions depend on assumptions about the state of the early universe before the hot thermal 
phase, for instance the type of inflationary phase, and they are therefore strongly model 
dependent. 
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6 Outlook 



The main result of the paper is the production rate of gravitinos for supersymmetric 
QCD at high temperature to leading order in the gauge coupling. The result is valid for 
gravitino masses larger or smaller than the gluino mass. 

As expected the gravitino production rate depends logarithmically on the gluon plasma 
mass which regularizes an infrared divergence occuring in leading order. Following the 
procedure of Braaten and Yuan, the result is obtained by matching contributions to the 
gravitino self-energy with soft and hard internal gluon momenta and by using a resummed 
gluon propagator for the soft part. As a byproduct a new result for the axion production 
rate in a QED plasma is obtained which is slightly smaller than a previously published 
result. 

The QCD coupling is large, and even at temperatures T ~ 10 10 GeV the usually 
assumed separation of scales, g 2 T <C gT <C T appears problematic. Hence, higher-order 
corrections to the gravitino production rate may be sizeable. Further, it is of crucial 
importance to gain some understanding of the influence of the non-perturbative magnetic 
mass scale ji g on real-time processes in general. 

The thermal gravitino production rate plays a central role in cosmology since it is 
closely related to the dark matter problem. For many supersymmetric extensions of the 
standard model this rate defines a limiting temperature beyond which the standard hot 
big bang picture becomes inconsistent. At present supersymmetric theories offer several 
interesting candidates for cold or hot dark matter. It is an intriguing possibility that the 
gravitino itself is the dominant component of cold dark matter. 

We would like to thank T. Asaka, O. Bar, D. Bodeker, O. Philipsen and M. Plumacher 
for helpful discussions. The work of A.B. has been supported by a Heisenberg grant of 
the D.F.G. 
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Appendix A 



In the following we shall derive the prefactor of the self-energy, 

E(P) oc n 
extending the discussion in section 2. 



(a.i; 




(a) 






Figure 6: Contributions to the gravitino self energy. 



The gravitino self-energy takes the form (cf. eq. (|TT|) ). 



oc -tr 



n^(p),s(P)^...^(P) 



2ml 



3m%M 2 



tr 



[YnWP^iK,) . ..F Ta (K 1 )T(P 1 )[r,r] 



(A.2) 



The one- and two-loop contributions for the pure gauge theory in resummed perturbation 
theory are depicted in Fig. ([6]). The contributions (a) and (b) represent for the gluino 
line the first two terms of the gluino propagator. This corresponds to the substitution, 



A(P X ) . . . A(Pi) -> PiA(P 1; v) +VP(P, v) . 



(A.3) 



This is the general form of the gluino propagator for m~ g = because of chiral symmetry 
and the fact that the velocity v appearing in the gluon propagator is the only other vector 
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available apart from the momentum P\. With 

-f-flu = -2Y (A.4) 

one then reads off the factor flA.li for the contributions (a) and (b). The same arguments 
apply for the contribution Fig. (6d). 

For Fig. (6c) one obtains for the gluino line, 

A(P 1 )...A(P 1 )^f l7 T 2 - (A.5) 

With 

7f i7T27m = -2f 2 7Ti , (A.6) 

the interchange Pi <-> P2 and the property of the gluon propagator A(— K) = A(K) one 
obtains the factor (1A.1I) also for this contribution. 



Appendix B 



In this appendix we explain the calculation of the contribution of hard virtual photons 
to the axion production rate T^ a,Td (E). We start by reconsidering the defining equation 

m 



n B (E)T^ d (E) 



2 , d % 1 



4vr 2E 



d 3 Pi 



S (2vr) 3 2P,J 



(2tt) 4 5 4 (Pi + P 2 - P - P 3 



nl?I\M\ 2 e(\ Pl -p 3 \-k cut ) 
where the matrix element squared for e _ 7 — > e~a is given in Eq. 



and 



n 



FBF 
total 



n F (E 1 )n B (E 2 )(l-n F (E 3 )). 



(b.i; 



(B.2) 



A great simplification is achieved in the computation of (IB. II) if one uses as reference 
momentum the difference vector k = pi — P3, i.e., we write 

d 3 Pl 



2E X 



8(P^)e(Ei)dE 1 d 3 p 1 

J d 3 k5 3 (k + p 3 - Pl )5(P 1 2 )0(P 1 )rfP 1 d 3 p 1 
S(E 2 - |k + p 3 \ 2 )Q{E 1 )dE 1 d 3 k. 



(B.3) 



Further, 

d 3 p 2 
2Eo 



S\P 1 + P 2 -P-P 3 



5(P 2 )Q(E 2 )d i P 2 5 4 (P 1 + P 2 - P - P 3 ) 
5((P + P 3 - P x ) 2 - (p - k) 2 )0(P + P 3 - Pi).(B.4) 
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We now use rotational invariance to choose 



k = k (0,0,1), 

p = P (0, sin 9, cos 9), 

p 3 = P 3 (cos^sin^, sin^sin^, cos 9), (B.5) 



which implies 



and 



s = (P 1 + P 2 ) 2 = (P + P 3 ) 2 = 2PP 3 (1 -sin 9 sin <p sin 9- cos 9 cos 0), 

t = (Pi - P 3 ) 2 = (£?!- - fc 2 , (B.6) 



|k + p 3 | 2 = El + A; 2 + 2P 3 A; cos 9, 

|p-k| 2 = E 2 + k 2 -2Ekcos9. (B.7) 



It follows that 



2£;P \ 2kE 



5(P 2 -|k + p 3 | 2 ) = J— 6 (cos 6- *% 5 (B.8) 
v 1 1 ^ ' 2kE 3 V 2£P 3 J K ' 

The integrations over the (^-functions yield the following 0-functions (where we use also 
the 6-functions 6(Pi), 6(P + E 3 - P x ) and E = p > 0, P 3 = p 3 > 0: 

1. ) From the integration over cos 9 we get: 

cos# < 1 — > k > Pi - P 3 , 
cos#>-l -> Pi>|P 3 -A;|. (B.9) 

The second of these constraints is equivalent to 

P 3 - Pi < k < Pi + P 3 . (B.10) 

2. ) From the integration over cos# we get: 

cos#<l -> |P- fc| < P + P 3 -Pi, 
cos#>-l -> k>E 3 -E 1 . (B.ll) 

The first of these constraints is equivalent to 

Pi - P 3 < k < 2P + P 3 - Pi. (B.12) 
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After integrating out the ^-functions we therefore have: 

1 1 



n B {E)T^\E) 



(B.13) 



2 7 7T 4 E 2 

where Q is the product of all 0-functions that restrict the integrations over Ei, E 3 and k, 

n = e(k-k cnt )e(k-\Ei-E 3 \) 

Q(Ex + E 3 - k)G(2E + E 3 - E 1 - k) 

e(E 1 )B(E 3 )e(E + E 3 - Ex) . (B.14) 
Since only s depends on 0, we can integrate out also this angle without difficulty: 



#|M| 2 = j 2 



-2s 2 



2s -t 



3e 2 vr 



(Ei — E-. 



) 2 -k 2 



2f 2 

(E 3 + E 1 ) 2 (E + E 2 ) 2 
k 4 



-1 + 



= 9- 



2 E 2 + E 2 + 2EE 2 



(B.15) 



We thereby obtain the result given in Eq. (1271) in the main text. We now rewrite the 
expression for ( IB. 14}) using 

e(E 1 +E 3 -k)=l-e(k-E 1 -E 3 ). (B.16) 

We use e(Jfe -Ei- E 3 )Q{k - \E X - E 3 \) = G(k -Ex- E 3 ) and thus get 

Q = \e(k - k cut )Q(k - \Ei - E 3 \)e(2E + E 3 - Ei - k) 
- &(k - k cut )Q(k -Ei- E 3 )Q(2E + E 3 - E x - k) 
x Q(E 1 )Q(E 3 )Q(E + E 3 -E 1 ). 

We multiply the second term in the brackets of Eq. (1B.17j) with 1: 

1 = 6(A; cut -Ei- E 3 ) + Q(Ei + E 3 - k cut ), 



(B.17) 



(B.18) 



and note that 



Q{k - k cut )Q(k -Ei- E 3 )Q(k c 

eik-k^eik^-Ex-Es), 



Ei-E 3 



(B.19) 



and 



G(k - k cnt )Q(k -Ei- E 3 )Q(Ei + E 3 - k cut ) 
Q(k-Ei-E 3 )Q(E 1 +E 3 -k cvlt ). 
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^B.20) 



The contribution from the first term on the r.h.s. of Eq. (1B.18j) is zero in the limit 



^cut — ■ ► 0. We see this by integrating over k from k cnt to 2E + E 3 — E\. The resulting 
expression has terms ~ lALt? ~ V^cut- Since from k cut > E% + E 3 it follows that both 
Ei and E 3 are smaller than k cnt it is easy to see by power counting that the expression 
after the k integration is of order k cnt . Then we are left to consider: 

n B (E)T h ™ d (E)=g l + g 2 , (B.21) 

where 

9i = y^&l dE ^J dE 1 n™£G(E + E 3 -E 1 ) 

x J dk&(k - k cnt )Q(k - \E X - E 3 \)Q(2E + E 3 -E 1 ~ k)g, 

WE* Jo dEs I dElU ™£ &i - E + Es ~ E ^ & ^ + E *- k ^ 
J dk®(k -E x - E 3 )Q(2E + E 3 - E x - k)g, (B.22) 

The integral over k in g 2 is nonzero only if 

E x + E 3 < 2E + E 3 - E x E x < E. (B.23) 

In the limit k cnt — > we therefore get: 

92 = Y^p^ J, dE * J, dE ^<otI ( E i - E ) l( E + E i) E z + E l {E- E^] , (B.24) 
We rewrite this result for later use as follows: 

92 = 77 n J9 T^r / dE 3 / dEx In 



92 ~ 2¥ 



x 



167r 3 f 2 E 2 Jo Jo ' \ E 3 

d 
~dE~x 



x Q(E - SO^r \n^{ElEl - E 2 E^) 



(B.25) 



We now turn towards the computation of g\. We first multiply by 1: 

9i = 9MK ut - \Ex - E 3 \) + ^6(1^! - E 3 \ - k cut ) = gil + g 12 (B.26) 



Note that 



G(k - A; cut )6(A; - \Ex - E 3 \)Q(k cut - \Ex - E 3 \) 

G(k - k cnt )Q(k cut - \Ex - E 3 \) (B.27) 
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and 



G(fc - A; cut )6(A; - \E, - £ 3 |)©(|£i - E 3 \ - fc cut ) 
= Q(k - \E, - E 3 \)Q(\E! - E 3 \ - k cnt ). (B.28) 



We therefore have 

1 I 



i i roo poo 

911 = 1 dEs Jo dElH ^ Q{E + E3 ~ E ^ k ^ ~ ^ " E ^ 

x J dkQ(k - k cnt )Q(2E + E 3 - E l - k)g (B.29) 



1 1 r°° 

9l2 = 



2 7 7T 4 E 2 



rOO pOO nnn 

/ dE 3 / dE in ™£ Q(E + E 3 - E 1 )Q(\E 1 - E 3 \ - k cut ) 
Jo Jo 



X 



J dkB(k - \Ex - E 3 \)Q(2E + E 3 -E l - k)g. (B.30) 



Consider first g n . The integration of g over k can be carried out easily. We do not write 
down the result explicitly but note that it contains terms ~ V^cut an d ~ 1/^cut- The 
integration over E 1 is done next. From A; cut > \E 1 — E 3 \ we get 

E 3 - k cnt <E 1 <E 3 + k cnt . (B.31) 

In the limit k cut — > we can therefore set E\ = E 3 in the distribution functions n^f and 
get 

e 2 

9u = 



poo 

in B (E)J dE 3 E 2 n F (E 3 )(l-n F (E 3 )) 



3tt 3 / 2 

= ^n B (E). (B.32) 
Now we turn towards the computation of gu- We insert 

1 = Q(E 1 -E 3 ) + Q(E 3 -E 1 ). (B.33) 
Then g 12 = g 12 i + ^122 with 

i 1 poo poo 

9121 = 2vWo dE3 i dE ^ME + E 3 -E 1 )Q{E 1 -E 3 -k cnt ) 
x J dk&(k-E l + E 3 )&(2E + E 3 -E 1 -k)g, 

I i /-00 poo 

9122 = ^7-4^ J q dE 3 j o dE 1 n™£e(E 3 -E 1 -k cut ) 

x J dk&(k-E 3 + E 1 )Q(2E + E 3 -E 1 -k)g. (B.34) 
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The integration over k gives 



0121 
0122 



167T 3 / 2 E 2 Jo 

e 2 1 r°° 
~16tt 3 / 2 ^2 J 



dE 3 
dE 3 



E 3 +E 



«^l"<total 



E 3 +k cut E% — E 3 

E3-k cut FBF E 2 [E\ + E 2 ) 

"^"total 







£"i — E 3 



(B.35) 



The logarithmic dependence on k cut is extracted by a partial integration with f'(Ei) 
1/(E 1 - E 3 ), f(E 1 ) = In - £ 3 |/£ 3 ). 
For the surface term we get: 



g2 ™ 
^surface = - —^ n B {E) dE 3 \li 



'fceut\ E$eX V (E 3 /T) 

v E 3 J (exp(E 3 /T) + If 



lrp3 



24vrf 



In 



2T 



cut 



3 C ; (2) 
2 7+ C(2) 



;b.36) 



In the remaining term, which is given by — / dEif(Ei)g'(Ei), k cut can be set to zero. 
Writing 5-12 = fiWace + partial we obtain: 



fl'partial 



16vr 3 / 2 E 2 

e(s 3 -Si) 



d£ 3 / dEtBiE + E 3 - Ex) In 



l-^i — -^3 
E 3 



. d 




"Ue x L 


d 


n 


dE x 





„FBF p2f Tp2 , 77^ 
"total h 2 (A + h 3 



<°lE\El + El) 



(B.37) 



Performing the differentiation and combining the results for g\ and g2 leads to the final 
result Eq. ( 1291) given in the main text. 



Appendix C 



In this appendix we describe in some detail the calculation of the hard virtual gluon 
contribution and of the other non-singular contributions to the gravitino production rate 
phard^g^ ^y e s ^ aT ^ considering the defining equation (HI]) . By summing the cor- 
responding squared matrix elements of table 1 with the appropiate multiplicities and 
statistical factors, we get 



IM, 



BBF 



1 



3m l, 



2g 2 (N 2 



V 



M 2 



s + 2t + 



It 2 



(N + rif) + 2sn } 



(c.i; 
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IM, 



BFB 



1 + 



IM, 



m 9 

3m i 



First we note that since s 



4g 2 (iV 2 - 1] 
M 2 

4g 2 (iV 2 - 1 
M 2 

!-tiwe may write s 



2s 2 \ 

-t-2s- — UN + n f ) 



2tn, 




t 



— (N + Tlj 



(C.2) 
.(C.3) 



2t = t — u and 



t s + t 



t 



s 
u 



(C.4) 



The difference t — u and 1/t — 1/u is odd under exchanging Pi and Pi. If the remaining 
integrand and the measure is even under this transformation, the integral over such terms 
will be zero. Therefore in \Mbbf\ 2 , the contribution of s + 2t will give zero. Further we 
may trade s in \M BBF \ 2 with — 2t. In \M FFF \ 2 , we may likewise substitute 



_s* 2s 2 

t + s + t^ t 

Therefore only the following squared matrix elements have to be considered: 

2s 2 



(C.5) 



|Mi| 
\M 2 \ 

\M 3 \ 



-t - 2s 



t 



t. 



(C.6) 



and we replace the matrix elements in ( HIT ) by 



IM, 



BBF 



IM, 



BFB 



IM, 



FFF\ 



1 + 



1 + 



1 + 





Ag 2 (N 2 - 


1) 




M 2 




m? > 


4^ 2 (iV 2 - 


1) 


3m 2 , y 


M 2 




m? > 


4^ 2 (iV 2 - 


1) 


3ml J 


M 2 





iMal^A^ + n/) - 2\M 2 \ 2 n f 



\Mi\ 2 (N + n f ) -2\M 2 \ 2 n f 



Mi] 



IM, 



(N + n/ ) . 



(C.7) 
(C.8) 
(C.9) 



|Mi| 2 is the axion matrix element which has been discussed in appendix B. The different 
statistical factors in the case of gravitino production do not change the structure of the 
contribution from \Mi\ 2 as compared to the axion case. Again we can extract the loga- 
rithmic dependence on the cutoff k cut by a partial integration. In the case of BFB, the 
surface term contains an additional term which depends on the energy of the gravitino, 
see Eq. ( 10. 14[) below. The other two matrix elements do not induce a logarithmic depen- 
dence on the cutoff k cnt , i-e. one can set fc cu t = to compute their contribution to the 
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gravitino production rate. The contribution from |M 2 | 2 = t can be obtained easily using 
the same methods as in the axion case, where now no partial integration is needed. We 
obtain 



P 

BBF(BFB) 



dVt p 1 
17 2E 



n 



d 3 pi 



96tt 3 



o 



t\ (2vr) 3 2E ; 

rE+E-i 

dE 3 / dEiriBBF(BFB) 



(2tt) 4 5 4 (P 1 + P 2 - P - P 3 )n BBF{BFB) \M 2 \ 



x { Q(E - E 1 



E — Ei - 
] ~E T 



2E 2 + (3E 3 -E 1 )(E + E 1 ) 



- e(E 1 -E 3 )^{2E-E 3 + E 1 ) 
+ 6(P 3 -E^ (-3P 3 + 3Pi -2E) 



(CIO) 



To compute the contribution from |M 3 ' 2 



= — it is convenient to choose different coordi- 

O I s 

nates to perform the angular integrations, namely 



q = p + p 3 = g(0,0, 1), 
p = E (0, sin 9, cos 9), 
p 2 = E% (cos sin 6 1 , sin sin 6* , cosi 

"'hard / 



(c.n; 



The calculation of this contribution to 1^ (E) then goes along similar lines as for the 
axion, i.e. the integration of the angular variables cos 9, cos 9 can be trivially performed 
using the 5— functions. This leads to several constraints for the integration over q, which 
can be performed without problems. The final result is rather compact: 



T e/ S 

BBF(FFF) 



dVtp 1 
17 2E 



3 d 3 pi 

tk (2^) 3 2^ 

E+E-i 



(2 7 r) 4 «5 4 (P 1 + P 2 - P - P 3 )n BBF{FFF) \M 3 



I roo rE+E-i 
327 Jo dE3 i dE ^BE(EEE) 



'''' 1 6(P 2 - E 3 )^J^ [^(P 3 - E 2 ) + E(E 3 + E 2 )} ). (C.12) 



v P + P 3 ' ~ v * " y E 2 
The full result for r^ ard (P) can be written as in Eq. (j%3"|) with 

I BBF = 32tt 3 (e E ' T + l) [(N + n } )li' B s F - 2n f I 



t 

BBF 



(C.13) 



I B FB=T\N + n f ) 



Li 2 (-e" 



-E/T\ 



7T 



l + 81n(2)) 



64a* n/ (e E / T + l)l BFB ,(C.U) 
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W = -32vr 3 (e E ' T + l) (N + n/)/?^. (C.15) 

The hard contribution to the collision term can be obtained by a numerical integration. 
Adding all the contributions we finally find: 

Cf rd (T) = j ' -0^n F (E)T^ d (E) 

( m|_\ SCWjN 2 - 1)T 6 
^ 3ml) 32tt 3 M 2 

+ 0.5781ri/|, (C.16) 

which yields, after adding the soft contribution, our final result for the gravitino collision 
term (H4"j) . 



x \(N + n f ) 



In 



' rp2 ' 
, ^CUt , 



1.7014 
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Erratum 



As recently pointed out by Pradler and Steffen (Phys. Rev. D 75 (2007) 023509 



|hep-ph/0608344|), Eq. (C.14) has to be replaced by 

Ibfb = -\v 2 ln(2)T 3 (iV + n f ) - Q4n 3 n f (e E/T + l) I BFB . (E.l) 

The difference between Eqs. (C.14) and (E.l) is a non-vanishing surface term in the BFB 
case, which in the partial integration leading from (B.24) to (B.25) was overlooked. 

As a consequence, the numerical factor multiplying (N + n/) changes in Eq. (C.16) 
from 1.7014 to 1.8782. Correspondingly, also the numerical factor multiplying (N + nf) 
in Eq. (44), and the overall factors in Eqs. (47) and (48) change. The corrected equations 
read 



C 6 (T) = j ^n F {E)(Tf{E)+Tt\E) 



3C(3)^ 2 (iV 2 - 1)T 6 



3m% / 32tt 3 M 2 

(j * 



T 



2 ' 



hi I — I + 0. 1002 



ml 



(N + n f ) + 0.5781n/ \ , (E.2) 



Y — 14- 10-- ( Tr ) ( mGeV ) 2 2 (E 3 ) 

tt d h 2 = tUoYq (T)n ra( i (T) h 2 p~ x 

= ^(^^)(^m(M' . (E.4) 



10 10 GeV/ V m G J \ 1TeV , 

Compared to Eqs. (47) and (48), the gravitino abundance is increased by 30%. (Note 
that also the prefactor l/(2\/6) in Eq. (33) has to be replaced by l/(4-\/6).) 



We thank Josef Pradler and Frank Steffen for informing us about the error in our paper. 
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